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ABSTRACT 

The least squares collocation algorithm for estimating gravity 
anomalies from geodetic data is shown to be an application of 
the well known regression equations which provide the mean 
and covariance of a random vector (gravity anomalies) given 
a realization of a correlated random vector (geodetic data). 
It is also shown that the collocation solution for gravity anom- 
alies is equivalent to the conventional least-squares -Stokes' 
function solution when the conventional solution utilizes prop- 
erly weighted zero a priori estimates. The mathematical and 
physical assumptions underlying the least squares collocation 
estimator are described and its numerical properties are 
compared with the numerical properties of the conventional 
least squares estimator. 
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INTRODUCTION 


The problem of recovering mean gravity anomalies from geodetic data has re- 
ceived much attention. The conventional approach to the problem is to employ 
the discrete form of Stokes' function to obtain an equation of condition for a 
least squares estimate of mean gravity anomalies. If the best available a priori 
estimates of mean gravity anomalies are included with proper weighting in the 
least squares loss function, the resulting estimates are known to be optimal in a 
minimum variance sense. • This approach has been simulated for satellite per- 
turbation data by Hajela^ , for gradiometer data by Reed3 and by Argentiero and 
Garza-Robles^, and for altimeter data by Gopalapillai^ and by Argentiero (et al).^ 
For insitu data types such as altimetry and gradiometry it can be shown that the 
procedure can be used to estimate local blocks of anomalies in local blocks of 
data without significant aliasing. 

An alternative approach to geodetic data reduction problems called least squares 
collocation has been suggested by Moritz."^ It is claimed that the least squares 
collocation method is a more general and more powerful parameter estimation 
procedure than the classical least squares method It has also been as- 
serted that least squares collocation is the only parameter estimation method 
which permits the simultaneous and optimal processing of heterogeneous data 
types."- *2 Moritz" has recommended leastsquares collocation as the preferred 
method for estimating mean gravity anomalies from gradiometer data. Rapp"^ and 
Smith *•* studied the problem of applying the collocation technique to estimat- 
ing mean gravity anomalies from altimeter data. Rapp*^ asserts that be- 
cause of the deterministic nature of the Stokes' tunction, the use of conventional 
least squares techniques in estimating mean gravity anomalies from altimeter 
data can lead to false or misleading results . This criticism is repeated by Uotila.* ® 

In the following discussion the least squares collocation procedure is viewed as 
a direct application of the regression equations which provide the conditional 
mean and covariance of a random vector ^mean gravity anomalies) given a reali- 
zation of a correlated random vector (geodetic data). It will be seen that this 
approach permits both a fast and simple derivation of the collocation algorithms 
and a derivation of their precise relationship to the least-squares -Stokes' function 
algorithm for estimating mean gravity anomalies from geodetic data. Finally the 
mathematical and physical assumptions underlying the application of the least 
squares collocation algorithm to the estimation of mean gravity anomalies are 
discussed. 

THE LEAST SQUARES COLLOCATION SOLUTION 

In this presentation, least squares collocation will be treated as a parameter es- 
timation procedure in a finite dimensional real vector space. Since for all 
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practical applications a finite amount of data is processed to estimate a finite 
number of parameters, no useful generality is lost by this approach. 

Let {y'} be a set of geodetic observations. The problem is to obtain from such 
an observation set a "best" estimate of a set of mean gravity anomalies {6g}. 

We will define the "best" estimate to be the conditional expectation of {5g} given 
a realization of the observations { Y'} . Since the smallest second moment of a 
random variable is the second moment about the mean, this is equivalent to ap- 
plying a minimum variance criterion. 

The starting point of the least squares collocation approach to obtaining the best 
estimate of {fig} is the assumption that one has full knowledge of the second order 
statistics of the anomalous potential field everywhere on and outside the reference 
geoid. (The first order statistics of the anomalous potential field are assumed 
to be zero.) Let P(xi) and P(X 2 > be the anomalous potentials at points xj and xi 
on or outside the reference geoid. We assume the possession of a function 
K (x 1 , X 2 ) such that 


E (P(X,)P(X2)) =K (X,,X2) (1) 

The function K(xi ,X 2 ) is the so-called covariance function, and it is generally 
defined to be invariant under rotations. Hence the second order statistics of the 
anomalous potential field are assumed to be independent of location. Let {Y| 
be a vector which is determined by the anomalous potential field and which, after 
suitable corrections for systematic error sources in the measurement process, 
is directly observable. Also, let {fig} be a set of globally distributed gravity 
anomalies. Since both {Y} and {fig} are determined by the anomalous potential 
field all second order statistics relating to the two random vectors can be readily 
derived from the covariance function. Hence define 

a)E(YYT) = A, b)E(figYT) = B, c)E(5gfigT) = C (2) 

Computational algorithms for obtaining matrices A, B, and C from a covariance 
function are developed by Moritz‘S and by Tscheming and Rapp.** The actual 
observations {Y'} obtained from the instruments are, of course, corrupted by 
noise. Hence 

Y' = Y+v, E(v) = 0, E(vi;T) = q (3) 

Equations (2) and (3) permit us to write the joint covariance matrix of the random 
vectors {Y'} and {fig} as 


2 


cov 


(4) 


’fig" 


■ C 

B 

Y' 


BT 

A+Q 


A realization of the random vector {Y'} is obtained by means of the actual meas- 
urements. Symbolically we do not distinguish between this random vector and its 
realization. We desire the conditional expectation and the conditional covariance 
of {6g} given a realization of the correlated random vector {y'} . By .assuming 
either that the random vectors are normally distributed or that the conditional 
expectation of {6g} is a linear function of the measurements we can resort to the 
familiar regression equations for the conditional mean and conditional covariance 
of a random vector given a realization of a correlated random vector.* The re- 
sults are 

6g = B(A + Q)-‘ Y' (5) 

COV [6g = C - B (A + Q)-* BT (6) 

The solution represented by Equation (5) is the least squares collocation estimate 
of a global set of gravity anomalies given a covariance function .and given the 
measurement set {y'} . 


In actuality, one would not attempt to estimate a global set of anomalies from a 
set of geodetic observations obtained from a certain area. It is only possible to 
significantly improve knowledge of gravity anomalies in the area covered by the 
observations. Decompose {6g} as follows: 



(7) 


where {5gi} is the set of anomalies covering the region where the measurements 
are available and where {^g 2 } is the set of anomalies outside of this region. Then 
the matrix B can be decomposed 


B = 


B2 


( 8 ) 


where 


B, =E (6g,YT), B2 =E (6g2YT) 


The least squares collocation estimate for {8g|} becomes 


6gi =B, (A+Q)-' Y' 


(9) 


( 10 ) 
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THE CONVENTIONAL LEAST SQUARES SOLUTION EMPLOYING STOKES’ 
FORMULA 

The discrete form of Stokes' formula provides a linear relationship between the 
value of the anomalous potential field at any point and a set of globally distributed 
mean gravity anomalies. Since {y} is a vector which is determined the anom- 
alous potential field, Stokes' formula provides a relation between {Y} and {6g} 
which after suitable linearization can be written as 

Y = S6g (11) 

The elements of the matrix S are obtained by evaluating Stokes' function at the 
required computation points. 

Equation (11) can be used as an equation of condition for a least squares estimate 
of {5g) . But the resultant solution would not be optimal unless all information 
were used. Consequently if one accepts the validity of a covariance function it 
would be proper to utilize the zero vector as an a priori estimate of {5g} with 
a weight provided by the inverse of the covariance matrix of Equation (2) c). The 
resultant loss function to be minimized has the form 

L (6g) = (Y' - S 6g)T Q-' (Y' - S 6g) + 6gTC-> (12) 

The estimator which minimizes th right side of Equation (12) is 


6g = (STQ-Is + C->)''STQ-> Y' (13) 

Equation (13) provides the standard least squares solution for {6gJ using Equa- 
tion (11) as an equation of condition and using the zero vector weighted accord- 
ing to a covari£mce function ns an a pr.ori estimate. 


DERIVATION OF AN EQUIVALENCE RELATION 

It cim be shown that the conventional least squares estimate of I^g) as defined 
by Equation (13) and the least squares collocation estimate of |6gJ as defined by 
Equation (5) are equivalent. Equation (11) defines the zero expectation random 
vector {Yj in terms of the zero expectation random vector »6gJ. Thus the co- 
variance matrix of !YJ and the joint covariance of {Yj and (6gl cim be obtaine<J 
in terms of the covariance matrix of (6gJ . Equation (2) c) provides the covari- 
ance matrix of {6gl as derived from the covari.ance function of Equation (1). 
liquation (11) along with Equations (2) c) juid (3) permit us to write 


4 


cov 


(14) 




■ C 

CST 

Y' 


SC 

SCST + Q 


The regression equation can again be used to obtain the conditional expectation 
of 5g as; 

6§ = CS'r (SCS'T + Q)"‘ Y' (15) 

A comparison of Equation (4) with Equation (14) yields 

B = CS’T , A = SCST (16) 

Hence the estimate of {6g} provided by Equation (15) is equivalent to the least 
squares collocation estimate of Equation (5). We can use the well known Shure 
matrix identity to translate Equation (15) into the alternative form: 

= (STQ-«S + C-' )■' SiQ-' Y' (17) 

Equation (17) is identical to Equation (13). This demonstrates that a standard 
least squares approach to estimating gravity anomalies from geodetic data which 
utilizes an a priori estimate weighted according to a covariance function yields 
a solution identical to what is obtained through least squares collocation. 


THE DERIVATION OF COVARIANCE FUNCTIONS 

The implementation of the collocation technique depends explicitly on the existence 
of a covariance function which provides the second order statistics of the anom- 
alous potential field. Either of two basic methods may be employed in deriving 
such a covariance function. Both methods assume that the statistical behavior 
of the anomalous potential field is independent of location. These methods are 
summarized below; a detailed exposition may be found in Reference 15. 

One method derives the model from a set of ground -based gravitational anomalies, 
as nearly "global" in scope as is available. The covariance function for me:m 
free air gravity anomalies is estimated as 

C(^) = \ (18) 

La N 

where 5g(0,^) and fig(0',X') are pairs of mean gravity anomalies in e<iual area 
blocks separated by a spherical distance \}/ and N is the s”m of products at that 
(iistance. The covariance function for the anomalous p' . '-ntial c.on then lx? ob- 
tained as: 
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K (X,, X,) = 


(8 


+ 1 I — 

^ y 1*1^ (cos J C( ^(/')l*J;(COS i//') sin i//' d (19) 


where X| and x> are two points on or outside the reference geold which are sep- 
arutcnl by spherical distance C^i Pv Is the Legendre polynomial function Jind 0 
and s are in essence geometric factors. 


This approach implies that the data set is a random sample of the Karth's gravi- 
tational features. In particular, there is an implicit assumption that the char- 
acteristics of gravity anomalies over ocean areas are the same as on land, or 
that the sample contains a proportionate distribution of anomalies over hmd and 
ocean areas. 

An alternative approach is to obtain the covariance function for the anomalous 
potential directly from a model for the spherical harmonic coefficients of the 
potential field. Let o?(8) be a model for the second moment alxiut zero of nor- 
malized spherical harmonic coefficients of degree 8 . Then the covariance func- 
tion for the anomalous potential can be written as 


K(x,.x,)= s'*' (28+ l)oJ(q^, Sj^lP.costi;/) (20) 


0 


where again, s is a geometric factor, and [j/ is the spherical distance between 
points X| and X 2 . One model for the second moment of spherical bnrmonic 
coefficients is the well known "Kaula rule of thumb." 



Other models are described in Reference 15. 

The implicit assumption of this procedure is that the magnitude of spherical 
harmonic coefficients decreases monotonically with degree and is representable 
as a simple function of degree. There are indications that the real field may 
not have this property’. Also the covariance functions derived from the alter- 
native procedures outlined above can differ substantially.'^ 
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COMMENTS 


The least squares collocation algorithms are a form of the regression equations 
which provide the conditional mean and covariance of a random vector (mean 
gravity anomalies) given a realization of a correlated random vector (geodetic 
data). Tapleyi^ has showm that the collocation model cao be cast into a form in 
which standard least squares reduction procedures are applicable and that the 
solution so obtained is identical to the collocation solution. In addition, the 
conventional approach to estimating mean gravity rjiomalies from geodetic data 
which relies on Stokes' formula and a least squares estimation algorithm reduces 
to the collocation model provided an a priori estimate is weighted according to a 
covariance function and is included in the loss function. Hence, assertions that 
the least squares collocation procedure represents a more general and more 
powerful method for estimating mean gravity anomalies than what is provided by 
the conventional approach arc unfomided. 

Several investigators’*'^ have found that the results of a least squares colloca- 
tion are not sensitive to the choice of a covariance function. The results of this 
paper provide an explanation of this fact. la effect, the covariance function 
supplies a weight for a zero a priori estimate of mean gravity anomalies in a 
least squares reduction procedure. In any meaningful estimation process the 
resultant solution is dependent mainly on the information content of the data 
rather than the quality of the a priori estimate. Hence, in a practical simula- 
tion or real data reduction one should expect the results of the collocatlDn pro- 
cedure to be relatively insensitive to the choice of a covariance function. Fur- 
thermore, this implies that for most applications the two estimation procedures 
should yield almost identical results. 

If the a priori information provided by a covariance function is to be included in 
a solution for mean gravity anomalies it is better numerically to do so by means 
of the Stokes' formula and the conventional least squares method since this pro- 
cedure involves the inversion of a matrix whose dimension is the number of es- 
timated parameters. The collocation algorithm as given by Equation (10) implies 
the inversion of a matrix whose dimension is the size of the data set. 
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